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ABSTRACT 

We have mapped the key mid-IR diagnostics in eight major merger systems of the Toomre Sequence 
(NGC 4676, NGC 7592, NGC 6621, NGC 2623, NGC 6240, NGC 520, NGC 3921, and NGC 7252) 
using the Spitzer Infrared Spectrograph (IRS). With these maps, we explore the variation of the 
ionized-gas, PAH, and warm-gas (H2) properties across the sequence and within the galaxies. While 
the global PAH interband strength and ionized gas flux ratios ([Ne III]/[Ne II]) are similar to those 
of normal star forming galaxies, the distribution of the spatially resolved PAH and fine structure line 
flux ratios is significant different from one system to the other. Rather than a constant H2/PAH 
flux ratio, we find that the relation between the H2 and PAH fluxes is characterized by a power law 
with a roughly constant exponent (0.61 ± 0.05) over all merger components and spatial scales. While 
following the same power law on local scales, three galaxies have a factor of ten larger integrated (i.e. 
global) H2/PAH flux ratio than the rest of the sample, even larger than what it is in most nearby 
AGNs. These findings suggest a common dominant excitation mechanism for H2 emission over a 
large range of global H2/PAH flux ratios in major mergers. Early merger systems show a different 
distribution between the cold (CO J=l-0) and warm (H2) molecular gas component, which is likely 
due to the merger interaction. Strong evidence for buried star formation in the overlap region of the 
merging galaxies is found in two merger systems (NGC 6621 and NGC 7592) as seen in the PAH, [Ne 
II], [Ne III], and warm gas line emission, but with no apparent corresponding CO (J=l-0) emission. 
The minimum of the 11.3/7.7 /jm PAH interband strength ratio is typically located in the nuclei of 
galaxies, while the [Ne III/[Ne II] ratio increases with distance from the nucleus. Our findings also 
demonstrate that the variations of the physical conditions within a merger are much larger than any 
systematic trends along the Toomre Sequence. 

Subject headings: galaxies: evolution, galaxies: interactions, galaxies: starburst, ISM: molecules, 
galaxies: individual (NGC4676, NGC7592, NGC6621, NGC2623, NGC6240, 
NGC520, NGC3921, NGC7252 



1. INTRODUCTION 

Interactions and mergers are important drivers of 
galaxy evolution and are responsible for the formation 
of the most luminous infrared galaxies observed locally 
and many at high redshift. In particular, major merg- 
ers can completely disrupt the stellar and gaseous mor- 
phologies of the merging galaxies, transform spirals into 
massive ellipticals, and fuel both po werful starbursts and 
massive nuclear black holes (e.g., I Barnes fc Hernquistl 
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1992). In the most extreme cases, mergers can become 
Ultra-Luminous Infrared Galaxies (ULIRGs) which are 
in many cases powered by massive starbursts with star 
formation rates sufficient to form the entire stellar pop- 
ulation of an L* galaxy in only a few dynamical times 
(lArmus et all 119871: iSanders et al.lll988t faeckman et al.1 
119901: IMurphv et al 1119961) . While most mergers do not 
produce ULIRGs, they invariably trigger large-scale star- 
bursts and lead to a large-scale redistribution of the 
progenitors' gaseous and stellar disks. At high red- 
shift, mergers are responsible for the rapid evolution 
of the IR luminosity function seen in deep imaging 
surveys, and the formation of some of the most lu- 
minous galaxies d etected at both UV and far-infrared 
wavelengths (e.g.. iSteidel et al.| [l996t iBlain et "all 120021: 
iChapman et al.ll2005l ). At low redshifts. interactions and 
mergers provide us a detailed picture of the most rapid 
and intense star formation known, and the generation 
and fueling of active nuclei. 

The last decades have brought immense progress 
in the observations and modeling of both the stel- 
lar and the gaseous components of interacting galax- 
ies. F rom the pioneering s tudies of iToomre &: Toomrd 
(I1972T). and ISchweizerl (1 1 982f ) through the work 
of IBarnes fcH crnquist (1992) to the more recen t 
work of iFerrareseet all (120061). iCoteetall (l2007h , 
iKormendv et all (|2009f l and iHopkins et all (|2009[ ). it 
is now clear that interactions and mergers produce 
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large-scale (tens of kpc) stellar tidal tails, distorted 
and /or overlapping disks, central starbursts, and in 
some cases, stellar remnants that have the properties 
of elliptical galaxies. The gas, being dissipative, can 
rapidly form dense molecular disks , fueling both star 



bursts and an active nucl e us (e.g.. iScoville et al.1 



Mihos fc Hernquistl Il994t iCharmandaris et a! 



1991 



2002; 



Armus et al.l 120041 ). Hubble Space Telescope observa 



tions of nuclei in merging galaxies have shown an in- 
creasing nuclear luminosity density in the o ptical /near- 
infrared light with advancing merger stage (ILaine et al. | 
2001 iRossa et al.ll2007t iVeilleux et all l200l iHaan et al.l 
2011D . suggesting an increase of nuclear starburst activ- 



ity. However, the fate of the dust in galactic interactions 
and mergers, and its direct relation to the gas and the 
young stars is much less well understood. For instance, 
a change in dust grain size, ionization stage, neutral to 
ionized gas fraction, or gas temperatures as function of 
merger stage could be expected. The dust responsible for 
re-radiating much of the UV light into the far-infrared 
can optically obscure the most energetic regions of star 
formation in interacting galaxies, which makes it chal- 
lenging to explore the physical connection between the 
stars, gas and dust in the nuclei, disks, and overlap re- 
gions where the brightest, most intense starbursts are 
triggered. To address the physical conditions in major 
mergers and excitation mechanisms such as shocks and 
starbursts, it is necessary to resolve the spatial varia- 
tions of the warm dust, molecular gas, and ionized gas 
on subgalactic scales. 

Here we present Spitzcr IRS spectral mapping observa- 
tions of a sample of eight nearby, IR-bright galaxies that 
span the range from early through mid- to late-stage ma- 
jor mergers. Our maps of the mid-IR features reveal the 
spatial variations of gas temperatures, masses, dust grain 
sizes, and the ionization state of the interstellar medium. 
In § [5] we describe our sample, the observational setup, 
and the data reduction. The extracted spectra, mid- 
IR feature maps, dust and warm gas properties (e.g., H2 
temperature maps and masses) are presented in §E1 com- 
pared to normal star-forming galaxies, and discussed in 
the context of the physical properties of the ISM. Finally, 
we provide a brief overview of our conclusions in § 21 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Sample 

By e mploying the Spitz er IRS in Spectral Mapping 
mode (|Houck et al.l 120041) . we measure the gas and 
dust properties within a set of eight interacting/merging 
galaxies representing the early (NGC 4676, NGC 7592, 
NGC 6621 - distinguishable parent disks and large tidal 
tails), middle (NGC 2623, NGC 6240, NGC 520 - 
highly overlapping disk and/or double nuclei) and late 
(NGC 3921, NGC 7252 - single nucleus, elliptical-like 
light distribution) stages of the major merger process. 
In total, we obtained 44.9 hrs of Spitzer observing time 
in cycle 2 (P.I. Lee Armus) to map these eight galaxies 
with the SL and LL modules of the IRS. All eight galaxies 
have Spitzer Infrared Array Camera (IRAC) and Multi- 
band Imaging Photometer (MIPS) imaging, and are, ex- 
cept for NGC 6240, part of t he classical "Toomre se- 
quence" of m erging galaxies ()Toomre fc Toomre] Il972t 
iToomrd ll977l) . Note that the classical "Toomre Se- 



quence" lists galaxies only in rough order of complete- 
ness of the merger, and is not meant to represent a pre- 
cise order in time sequence of the merging process. An 
overview of the sample is given in Table [TJ The galaxies 
are all nearby (mean distance of 79 Mpc) , IR bright (log 
Lir\Lq] =10.3-11. 971.). and have a wealth of available 
ancillary data, from radio through x-rays (HI, CO, HST 
WFPC2 and ACS, NICMOS, VLA, Chandra). To illus- 
trate the stellar distribution of these major mergers, we 
show in Fig. [JJ the IRAC 3.6 /im images of the sample. 

2.2. IRS Observations and Data Reduction 

The Short-Low (SL) and Long-Low (LL) modules of 
the IRS are used to create 3D spectral data cubes over 
the central 1-2 square arcminutes in each system (corre- 
sponding to physical sizes of ~8-28 x 22-110 kpc) and 
over the entire wavelength interval (5-38 /im). In all 
cases the maps are large enough to cover the entire main 
bodies (as defined by the optical and NIR images) and 
overlap regions. All SL and LL maps have half slit- width 
steps in the dispersion direction to avoid gaps and pro- 
duce filled maps. The SL integrations are all 60 sec with 
two cycles per position. The LL integrations are 120 sec 
(one cycle), or 30 sec (3-4 cycles) per position to max- 
imize Signal/Noise (S/N) and avoid saturation near the 
nuclei. The SL maps use sub-slit steps in the spatial di- 
rection to cover the intended area. We have designed the 
maps to cover the central regions of each galaxy, and at 
the same time provide multiple sky positions (from the 
non-primary slit) for background subtraction (see Fig. [TJ . 
The IRS data reduction was performe d using the interac- 
tive IDL software package CUBISM (jSmith et al.l 12007ft 
that takes IRS Basic Calibrated Data (BCD) produced 
via spectral mapping mode, and constructs spectral data 
cubes suitable for detailed spatial and spectral analysis. 
CUBISM allows the user to construct 2D line and/or 
continuum images, and also extract ID spectra at any 
spatial position in the map (see Fig. [lj. Background 
spectra are constructed from our sky observations off 
the galaxies (when the primary slit is mapping the tar- 
get, the secondary slit is looking at the sky next to the 
galaxy). The background is subtracted from the BCDs 
before the cubes are constructed. At a first step, line 
([Nell], [Nelll], H 2 (S(0)-S(7)), [Ar II], [S ill]) and con- 
tinuum (dust and PAH feature) maps and ID spectra are 
extracted and fit within CUBISM. These maps are later 
used as a benchmark for our maps created with PAHFIT 
(see § r~ 



11 This refers to estimates of the 8—1000 (im infra red luminosity 
constr ucted from IRAS flux densities as defined by I Sanders et all 
(2003) 
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3. RESULTS AND DISCUSSION 

To investigate the properties of the warm gas and dust 
in our galaxies, we have studied the spatial distribu- 
tion and global properties of the main mid-IR features, 
which are: the primary PAH complex features at 6.2 /im, 
7.7 fim, 11.3 jLtm, 12.7 /im, 16.6 /im, and 17 /im, the 
primary diagnostics for the ionization state of the low- 
density gas ([Ne II] 12.81 /im, [Ne III] 15.55 /im, [Ar II] 
6.99 /im, [Si II] 34.82 /im, [O IV] 25.9 /urn, and [Fe II] 
25.99 /im emission line J^l) . as well as the pure rotational 
lines of H2 S(0) to S(7) which probe the warm molecular 
gas (a detailed study of the H2 temperature and masses 
is elaborated in § I3.3.2j) . Two different approaches are 
conducted. First, we study the global properties of the 
merger components, which has the advantage of high 
signal-to-noise ratios and allows us to compare our re- 
sults with spatially integrated mid-IR properties of nor- 
mal star-forming galaxies. Second, to investigate in de- 
tail the distribution of the warm gas and dust, we map all 
mid-IR line features using a customized PAHFIT routine 
and also compare these maps with the cold gas distri- 
bution (using radio interferometer CO and HI emission 
lines) of these major mergers. 

3.1. ID Spectra and Spectral Features 

To extract high signal-to-noise spectra and line fea- 
tures, we have integrated the flux over those regions of 
each merger system that show significant emission in the 
spectral maps. These regions for each system (ranging 
from 2-20 kpc) are outlined on the IRS continuum maps 
in Fig. [2]. The continuum maps are constructed using the 
mean flux at 5.3 — 5.7 /im and 13.6 — 13.9 /im, which are 
wavelength ranges with no significant emission line fea- 
tures. Spectra were extracted from the IRS SL (5-15 /im) 
and LL (15-37 /im) modules in spatially matched extrac- 
tion regions using CUBISM . Each spectrum is composed 
of data from the four low-resolution orders of IRS: SL2 
(5.25-7.6 /im), SL1 (7.5-14.5 /im), LL2 (14.5-20.75 /im), 
and LL1 (20.5-38.5 /im). Slight mismatches between or- 
der segments, resulting from small residual photometric 
and astrometric uncertainties, were addressed by fitting 
the continuum near the overlap region, then scaling SL2 
to match SL1, LL1 to match LL2, and finally SL to match 
LL. The scaled spectra were then concatenated and aver- 
aged in the overlap region. Typical scaling adjustments 
were ~10%. To measure the line flux of the mid-IR fea- 
tures, we fit the co mbined SL/LL spectra using PAHFIT 
(jSmith et al.ll2007l) . This spectral fitting routine decom- 
poses IRS spectra into broad PAH features, unresolved 
line emission, and grain continuum, with the main ad- 
vantage to recover the full line flux of any blended fea- 
tures. Reliable redshifts of our objects are taken from 
NED. The best-fit solution of the observed spectrum is 
computed using a x 2 minimization method. The code re- 
turns the best-fit parameters for each of the components, 
including the line or feature fluxes and the optical depth 
of the silicate absorption feature at 9.7 /im. 

PAHFIT provides several options to take dust ex- 
tinction into account. We employed PAHFIT assum- 
ing a uniform foreground screen extinction (SCREEN 

12 Note that the [O IV] 25.9 /im and [Fe II] 25.99 /im emission 
lines are blended in the LL data 



model, i.e. a screen of dust between the galaxy and ob- 
server), dereddening the emitted line intensities with the 
extinction corrections estimated fro m the depth of th e 
9.7 /im silicate absorption feature (Smith et al. 2007). 
The choice of a uniform screen model (over, e.g., a fully 
mixed dust geometry model or a galactic center extinc- 
tion) is based on observations which have shown that 
a simple screen attenuation is able to account reason- 
ably well for the diffuse IS M attenuation in normal and 
starburst galaxies (e.g. iCharlot "fe Fall 200H iSmith et al.l 
120071 iWild et aLl[20ilb . The silicate extinction as given 
by the fitted t at 9.7 /im ranges from 0.08 to 1.4 
(mean Tg.7 = 0.7 ± 0.5, median: 0.6) for our sample, 
which is slight ly larger than in normal galaxies (typi- 
cally Tg.7 < 0.2 ISmith et al.ll2007) and starburst galaxies 



(median T9.7 = 0.23 lBrandl et al 



2006 ), but smaller than 



in LIRGs (median T9.7 = 1 .4 IStierwalt et al.l 120 111 ) and 
ULIRGs (median t 9 . 7 = 1.6 lArmus et al ll2007h . To esti- 
mate the effect of dust extinction, we fixed in test-runs 
the extinction optical depth to zero and compared both 
results (extinction and extinction-free) using the H2 line 
flux at 9.66 /im. Since the H2 line at 9.66 /tm sits on 
top of the silicate absorption feature, we expect that this 
line is most affected by the extinction correction. We find 
that the models with no extinction have ~20% lower H2 
fluxes than the SCREEN models, indicating fairly small 
silicate absorption. Only for NGC 6240 and NGC 520 
we find a three times smaller H2 flux with no-extinction 
than for the screen model, which is due to the large sili- 
cate absorption in these galaxies (T9.7 ~ 1.3). For a more 
thorough discussion of how PAHFIT determines quanti- 
ties such as the uncertainty in the integrated line flux, 
line FWHM, line equivalent width, an d the used mid-IR 
extinction law see ISmith et all (pOOl . 

The combined SL/LL spectra and PAHFIT model arc 
shown in Fig. [3] for each region. The line flux and errors 
of the main mid-IR features are listed in Table [5J The 
5-38 /tm wavelength range contains many important di- 
agnostic lines, such as unresolved atomic fine-structure 
lines of Ar, Ne, O, Si, Fe, covering a large range in ion- 
ization potential. In addition to the fine-structure lines, 
numerous PAH emission features are easily detectable in 
our spectra; the fluxes of the strongest features at 6.2, 
7.7, 11.3, and the 17 /tm complex are listed in Table [H 
We see a large range in global line intensities and PAH 
interband strengths for our galaxies as described in more 
detail in the following sections. Furthermore, the pure 
rotational lines of H2 [S(0)-S(7)] are often very strong 
in starburst galaxies (Table [3]) and are used in § 13.3.21 
to determine the temperature and mass of the (warm) 
molecular gas as listed in Table [4] 

3.2. Mapping Spectral Lines and PAH Emission 

In Fig. [4] we present maps of the main mid-IR fea- 
tures for the eight major mergers in our sample. These 
continuum-subtracted maps were created using our own 
customized code with PAHFIT as core-algorithm. At 
first, we scale SL2 to match SL1, LL1 to match LL2, and 
finally SL to match LL by fitting the continuum near the 
overlap region (similar to the ID spectra fitting proce- 
dure). The combined LL cube is rotated and interpolated 
in order to match the SL cube, conserving the flux. The 
scale factors are obtained using the entire range of pixels 
that match the LL and SL aperture, assuming a uni- 
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form scale-factor rather than independent scale-factors 
for each pixel. The aligned SL and LL cubes are con- 
catenated on a common wavelength grid and averaged in 
the overlap wavelength regions. The maps in the com- 
bined SL/LL cubes (5.3-38.5 /zm) are smoothed by a 3 
pixel x 3 pixel box (5.55" x 5.55"), conserving the flux, 
to increase the signal-to-noise ratio of the spectra. Then, 
for each pixel, the spectrum is extracted and PAHFIT is 
used to decompose it. Finally, we obtain maps of the 
properties of the fitted mid-IR features (i.e., integrated 
line flux, line FWHM, line equivalent-width, the uncer- 
tainty in the line flux, the fit to the entire spectrum and 
the fit to the continuum). 

As a consistency check, we compared our line flux maps 
of the strongest line features with continuum-subtracted 
line flux maps created with CUBISM (see § 12^2"]) , The 
results show a very similar distribution of the strong line 
features, but CUBISM is less sensitive to reveal faint 
extended emission in the outer disk. Overall, the main 
advantage of fitting each pixel is that faint local varia- 
tions can be taken into account that would not be visi- 
ble by subtracting only the continuum averaged over all 
the pixels. Another important advantage of this method 
over CUBISM is to fit the full features before creating 
the maps, which would otherwise tend to lose wings of 
the broad PAH features and having difficulties to de- 
blend the PAH features from each other, and other con- 
tiunuum features (e.g., the silicate absorption). In some 
mid-IR feature maps (Fig. [3J we find faint rings which 
are artifacts of the continuum subtraction (also present 
in CUBISM derived maps, although in a smaller degree). 

3.2.1. Continuum Fluxes 

To characterize the basic properties of the spectral 
continuum we have mapped the flux densities for three 
rest-frame wavelength ranges: 5.3-5.7, 14.75-15.25, and 
29.5-30.5 (im. These wavelengths were originally chosen 
to cover a large baseline in wavelength while being least 
affected by PAH emission features, silicate absorption, 
or strong emission lines. The flux densities derived for 
5.5, 15, and 30 /im are the median values in the above 
wavelength ranges, respectively. Although the contin- 
uum fit derived with PAHFIT is not affected by features, 
we have chosen these wavelength ranges to be consistent 
with CUBISM-derived maps and previous IRS studies for 
comparison. All of our galaxies show a decrease of the 30 
to 15 /im continuum flux ratio towards their center. In 
general, also a smaller 15/5 /im ratio in the galaxies' cen- 
ter is found, but in some cases we find strong spatial de- 
viations from the azimuthal averaged profile (NGC 2623, 
NGC 4676, NGC 520, and NGC 3921). 

The spatially integrated 30/15 /im flux ratio of our 
sample ranges from 5 (NGC 6621, NGC 3921, NGC 7252) 
to 16 (NGC 2623) with a mean value of 7.3, which 
is similar to measurements of starburst galaxies (mean 
30/15 /mi flux ratio of 8.56; iBrandl et all [20061 ) . The 
same behaviour is seen for the spatially integrated 
15/5.5 /im flux ratio which ranges from 1.5 (NGC 6621, 
NGC 3921) to 6.3 (NGC 262 3) with a mean valu e of 3.4 
(3.8 for starburst galaxies; IBrandl et al.l I2006I ). Two 
merger components in our sample show significantly 
smaller ratios than the ratios averaged over the en- 
tire merger system: NGC 4676 South (30/15 /im=2.2, 
15/5.5 /xm=0.8) and NGC 6621 South (30/15 /mi=2.0, 



15/5.5 ixm=0.5). 

3.2.2. Polycyclic Aromatic Hydrocarbons 

PAHs are thought to be responsible for a series of 
broad emission features that domi nate the mid-infrare d 
spectra of starbursts (see, e.g., IPeeters et al.l 120041 ). 
They are observed in a diverse range of sources, with 
their strongest emission originating in photo-dissociation 
regions (PDRs), the interfaces between H II regions 
and molecular clouds. The relative strength of the 
different PAH bands is expected to vary with the 
size and the ionization state of the PAH molecule 
(|Draine fe Lil I200H ISmith et al.1 [200l . The 6.2 /im 
and 7.7 /im PAH bands dominate the emission of ion- 
ized PAH, whereas the 11.3 /im PAH emission is more 
intense for neutral PAHs. Variations of PAH fea- 
tur es in nearby g alaxie s have been studied r ecentl y 
bv IGalliano et al.l (l200ll HI onso- Herrero et all (120101) 
IPereira-Santaella et al.l (|2010fl . and iDi'az-Santos et al.l 
(1201 ID . including the effects of the radiation field hard- 
ness and intensity, the size distribution of the PAH car- 
riers, extinction effects, and the dependency on the age 
of the stellar population. 

Our spectra in Fig. [3] show that the spectral shape 
is usually dominated by strong emission features from 
PAHs. Here, we concentrate our analysis on several of 
the strongest PAH features at 6.2 /im, 7.7 /im (which 
is three blended sub-features centered at 7.42, 7.60, and 
7.85 /im), 11.3 /im (two blended sub-features at 11.23 and 
11.33 /im), and around 17 /im (two blended PAH sub- 
features centered at 16.4 and 17.1 /im). Other PAH fea- 
tures can be seen in the spectra (e.g., at 8.6 or 12.7 /im), 
but they have in general lower S/N or are blended with 
other features. The fluxes of the PAH bands at 6.2, 7.7, 
and 11.3 /im are all affected by dust extinction, but the 
11.3 /im is likely more affected by the 9.7 /im silicate fea- 
ture than the 6.2 /im feature. The 7.7 /xm PAH feature 
is affect ed by extinction to a similar level as t he 6.2 /im 
feature ()Draine '¥11 2001: Be irao et al.ll200l . 

In Fig. [5] and Fig. [6] we present the PAH 6.2/7.7 /mi 
and 7.7/11.3 /im ratio maps for the Toomre Sequence. 
For most of our galaxies the PAH feature at 11.3 /im is 
more extended than the PAH feature at 7.7 /im (factor 
of ~2, ranging from 1.3-3.3), equivalent to an increase of 
the PAH 7.7/11.3 /im ratio toward s the center, which is 
in agreement with previous st u dies (IGalliano et al.ll2008t 
IPereira-Santaella et al.l I2010t IDiaz-Santos et all 1201 1|) . 
This is expected, since the 7.7 /xm PAH feature is related 
to the ionized PAHs (which are more likely located in 
the central starburst region) rather than neutral PAHs, 
which are traced by PAH 11.3 /im and are located more 
likely in diffuse regions throughout the galaxies. How- 
ever, for three galaxies (NGC 4676 North, NGC 6240, 
and NGC 2623) we find the opposite behavior, namely, 
a decreasing PAH 7.7/11.3 /im ratio towards the center, 
which has bee n also found in the center of weak AGN 
host galaxies (ISmith et all 120071: lO'Dowd et al l 120091: 
IDiamond-Stanic fe Riekd 120101 ISales et al.l 120101 e.g.,) 
and suggested to be attributed to the selective destruc- 
tion of the 7.7 /im PAH molecule with the increasing 
hardness of the radiation field. This picture is supported 
by our observations since two of the galaxies mentioned 
above (NGC 6240 and NGC 2623) are characterized as 
low luminosity AGN hosts (see Table [T]), wheras none of 
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the rest of the sample seems to host an AGN. We find 
for most of our galaxies a decreasing PAH 6.2/7.7 fim 
ratio towards the center (typically a factor of 2). This is 
rather surprising since both bands, 6.2 /im and 7.7 /urn, 
are related to ionized PAHs, and h ave shown only very 
small variations in previous st udies ((Galliano et al.ll2008b 
IPereira-Santaella et all 120101 ). Since the PAH 6.2 and 
7.7 /xm ratio decreases towads larger grain sizes (larger 
number of PAH carbon atoms), an intuitive explanation 
is that the PAHs with larger grain sizes are more likely 
located in the outskirts rather than in the ce nter, which is 
in ag reement with models of PAHs in PDRs (jDraine & Li 
2001) . A possible correlation between PAH grain size and 
ionization state is discussed in the next section. 

3.2.3. PAH interband strengths and ionization 

To characterize the distribution and population of the 
PAH grain sizes and ionization states for each merger 
system, we present in Fig. [7] the variation of interband 
strength ratios of the PAH features and [Nelll] / [Nell] 
ratios (integrated over apertures of 3 x 3 pixels, or 
5.6" x 5.6") where significa nt emission i s pres ent, and 
compared to PAH models bv lDraine "&Lil(l200l . In gen- 
eral, neutral PAHs have large ratios of PAH 11.3/7.7 /urn, 
while PAH ions have much smaller values. For both 
neutrals and ions, there is a regular progression of PAH 
6.2/7.7 to smal ler values with inc reasing number of PAH 
carbon atoms (jDraine "&Lil[200l . 

We find a large variety of PAH distributions, rang- 
ing from ionized to neutral PAHs over a large range 
in grain sizes (from a few ten to thousands of car- 
bon atoms). For instance, NGC 4676 (primarily ion- 
ized PAHs) and NGC 520 (primarily neutral PAHs) 
show a very concentrated distribution in the PAH in- 
terband strength ratios, while, e.g., NGC 6621 (ionized 
and neutral PAHs) and NGC 3921 (primarily neutral 
PAHs, including the nucleus) exhibit more widespread 
PAH populations. We find no significant correlation be- 
tween PAH grain size (PAH 6.2/7.7 /xm) and ionization 
state (PAH 11.3/7.7 /im). Instead, we find a large va- 
riety of PAH grain sizes for a given charge state, e.g., 
for NGC 7592, NGC 6240, and NGC 7252. NGC 3921 
is mainly populated by neutral PAHs. For most of our 
galaxies, the central apertures (center of galaxies) have 
low PAH 11.3/7.7 ^m ratios, representing typically the 
most ionized PAH population for each galaxy, as one 
would expect from a central starburst or a low-luminosity 
active galactic nucleus (AGN). This finding is in agree- 
ment wit h other recent studies o f PAH populations in 
ULIRGs (|Dfaz-Santos et al.l f20lT| ) which show that the 
PAH 11.3 /urn is spatially more extended than the PAH 
7.7 fan feature. 

3.2.4. The [Ne IIIJ/[Ne II] line ratios: Probing the hardness 
of the radiation field 

In general the ratio of the 12.81 /jm [Ne II] and 
15.55 /im [Ne III] lines provides a diagnostic of ioniza- 
tion that is relatively insensitive to extinction, gas-phase 
abundances, and gas temperature (being low excitation 
lines in warm gas). Since the two ionization potentials 
of Ne are 21.6 eV [Ne II] and 41.0 eV [Ne III], the mea- 
sured [Ne III]/[Ne II] line flux ratio is sensitive to the 
average effective temperature of the ionizing stars, and 



hence might provide a measure of the hardness of the 
radiation field. 

The spectral maps of the [Ne III]/[Ne II] line ratios for 
our sample are shown in Fig. [5J For all our galaxies the 
nuclei have the smallest [Ne III]/[Ne II] ratios (see also 
right panel of Fig. [7]), while the circumnuclear regions 
tend to show (a factor of 2-4) larger ratios. Intuitively, 
this is in contrast to the expectation of a decrease of the 
hardening of the radiation field with distance to the cen- 
ter due to the small number of high energy photons avail- 
able at large distances to ionize Ne + to Ne 2+ . However, 
our result is in agreement with other recent observations: 
IPereira-Sa ntaella et al.l (|2010f) have found that the min- 
imum of the [Ne III] / [Ne II] ratio tends to be located at 
the nuclei for a samp l e of 15 LIRGs observed with Spitzer 
IRS; iThornlev et al.l ()2000D have shown that H II galax- 
ies have a nuclear [Ne III]/[Ne II] ratio of a factor ^2-3 
smaller than those measured in the extra-nuclear regions. 
A similar behavior of icreasing [Ne III]/[Ne II] ratios with 
increasing galactocentric dista nces has been obser ved in 
the nearby galaxies NGC 253 dDevost et al.l 120041), M8 2 
(|Beirao et al.l 1200^ . and NGC 891 (|Rand et al.ll200l . 
Galactic HII regions follow the same trend, with HII 
regions at increasing distances from the galactic nu- 
cleus showing larger [Ne III]/[Ne II] ratios (jGiveon et al.l 
2002). Our results confirm this trend, even if dust ex- 
tinction corrections are taken into account as a function 
of position in the galaxies (using the SCREEN model 
in PAHFIT). Metallicity-related effects are unlikely to 
be the dominant cause of the relative increase of [Ne 
III] strength outside the central regions of these galax- 
ies (jSniiders et al.l 120071: IPereira-Santaella et al.l 120101) . 
A more likely explanation may be that a fraction of the 
most massive stars are hidd en in ultra-compact HII re- 
gions (Ri gbv fc Riekd 120041 )) due to high dust densities 
in the nuclear regions, leading, e.g., to an decrease of the 
[Ne III] / [Ne II] ratio by a factor of ~2 if the density in- 
creases from 10 3 to 10 s cur 3 (jSniiders et al.ll2007t ). An- 
other explanation for this behaviour might be attributed 
to the variation in the age of the stellar populations, with 
a larger contribution of the old stellar population in the 
nuclei (e.g. due to stellar bulges), subsequently decreas- 
ing the average stellar age (and ionization) towards the 
center. 

3.2.5. Comparison of Global Mid-IR Properties to Normal 
Star Formation Galaxies 

The SINGS galaxies, being mostly isolated, star- 
forming systems with a large range in IR power, are 
a perfect sample against which to compare the mid- 
IR properties of our major mergers. In particular one 
might expect to find different PAH populations between 
these two samples due to, e.g., a different interplay 
between radiation field and the ISM. Since the PAH 
11.3 [im feature is thought to be produced by neutral 
PAHs, wh ereas the 7.7 feature arises primarily fr om 
PAH ions (jAllamandola et al.lll999HDraine &T11I2001D . a 
change in the neutral fraction should change the observed 
flux ratio. Harder radiation fields are also more efficient 
at destroying small PAH grains, an effect that may dom- 
inate the center of our galaxies. For all of our objects the 
range of the integrated (i.e., global) PAH strength ratios 
6.2/7.7 /im and 11.3/17 fim are in the range of [0.2-0.4] 
and [1-2.5], respectively. We find a range in the ratio 
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7.7/11.3 PAH bands of [1-6] and a range in [Ne III]/[Ne 
II] of [0.1-1]. The range and distribution of the ratios 
PAH 6.2/7.7 fxm, 11.3/17 fim, and [Ne III]/[Ne II] is very 
similar to interband s trength ratios found in the nuclei of 
22 starburst galaxies dBrandl et al. 2006) and for a large 
sample of SINGS galaxies (jSmith et al.ll2007l ). including 
normal star formation galaxies and low-luminosity AGN 
(Seyfert and LINER). This suggests that the physical 
properties in major mergers are, at least on global scales, 
similar to those in normal starburst galaxies, but can be 
very different on local scales within the individual merger 
systems as shown in the following sections. However, we 
do not see any global ratios of [Ne III]/[Ne II] > 2 in 
our sample, which is the case for five lower metallicity 
gala xies with H II-dom inated nuclei in the SINGS sam- 
ple (jSmith et al.l 120071 ) as well as in most blue compact 
dwarfs and in g alaxies with a strong AGN component 
IHao et all IpOOfl ). 

3.3. Properties and Distribution of the Warm Molecular 
Gas Component 
3.3.1. The Distribution of the Warm Molecular Gas 

Excitation of the warm molecular gas observed in the 
rotational emission lines of H2 (S(0)-S(7)) can com- 
monly arise from at least three different mechanisms: 
UV-excitation in PDRs surrounding or adjacent to the 
H II regions; shocks that collisionally excite the H2 
molecules; or hard X-ray photons capable of penetrating 
the molecular clouds a nd heating larg e columns of gas 
(jLepp fc McCravl 119831: iMalonev et aflll996t) . Observa- 
tionally separating these primary excitation mechanisms 
is not trivial in other galaxies: high UV radiation as 
well as high- velocity gas motions are both present in the 
nuclei of active galaxies and, to a lesser extent, in star- 
bursts. UV radiation can originate cither from the ex- 
isting stellar population, from current star-forming pro- 
cesses, or the AGN accretion disk. On the other hand, 
shock excitation can originate from the central engine, an 
outflow due to the large number of supernova remnants, 
or, finally from gas motions and cloud collisions due 
to galaxy-galaxy i nteractions (jShull fc Hollenbach1ll978l : 
iDraine et al.l H983). In addition, the low rotational tran- 
sitions arise from collisionally excited warm gas which 
complicates the distinction even further. However, if 
H2 emission is excited mainly by UV-radiation in PDRs, 
the fluxes of H2 lines and PAH features should correlate 
closely, since PAH emission features arise from the same 
mechanism in PDRs. 

In Fig. |H]we have mapped the emission from several H2 
rotational lines for each merger system. Our comparison 
between the flux distribution of the strongest H2 emission 
line (the extinction corrected 9.7 /jm H2 emission line) 
and PAH (using the combined main PAH features in the 
SL wavelength range at 6.1, 7.7, and 11.3 /Ltm) shows that 
(1) H2 has typically a larger extent than PAH with an 
average ratio of H2/PAH radius of 1.5 (ranging from 1.1 
to 1.9, including PSF size corrections), (2) some galax- 
ies show a significantly different H2/PAH flux ratio be- 
tween different merger components, most significantly in 
the early merger systems NGC 4676 and NGC 7592 (see 
§ 13. 5J) , and (3) at least one of our galaxies (NGC 3921) 
shows signs of a tail of H2 emission not associated with 
a strong PAH feature (see Fig. [4|. 



Five galaxies (NGC 7592, NGC 6621, NGC 6240, 
NGC 520, NGC 2623) show in Fig. Han offset of the peak 
H2 emission at 9.66 /im from the stellar nucleus, which 
is not seen in any PAH or other H2 features (e.g., S(2) 
and S(4)). One possible explanation for this offset is an 
underestimation of dust extinction (silicate absorption at 
9.7 /iin) which is strongest at the center of the galaxy, and 
hence suppresses the apparent H2 flux at 9.66 /im. How- 
ever, even without the 9.7 /mi H2 emission line, we find 
spatial deviations from a constant relation between the 
warm molecular gas (H2) and PAH (as described above), 
indicating that, at least in some merger galaxies, different 
physical conditions are present in the ISM (e.g., electron 
density, gas tem perature, or UV radi ation field intensity; 
for a review see I Galliano et al. 2008). The nature of the 
variation in the H2/PAH ratio will be discussed in more 
detail in § l3~3~4l 

3.3.2. Spatially Integrated Temperatures and Masses 

Temperatures and masses of the warm molecular hy- 
drogen (T<2000 K) are derived from the pure rotational 
H2 lines. Fig. [TUl shows the excitation diagrams using the 
extinction corrected line fluxes (see Table [3]) for different 
regions of the merger systems. In all cases an ortho- 
to-para ratio of 3 is assumed (Rigopou lou et al.1 120021 : 
iRoussel et al.l l2007) . The excitation diagram is basically 
a plot of the natural logarithm of the column density N$ 
divided by the statistical weight g{i) in the upper level 
of each transition against t he energy level Ej. A detailed 
description can be found in IRoussel et al.l (|2007[ ). includ- 
ing a discussion about the implications of using different 
ortho-to-para ratios. 

The excitation temperature of the line-emitting gas is 
the reciprocal of the slope of the excitation diagram. If 
the H2 gas is characterized by a single temperature (and 
assuming a single ortho-to-para ratio, the data will lie 
on a straight line, corresponding to a single component 
fit. We find that a simple single temperature fit does 
not work for most of our galaxies, since our data show 
a change of the slope with energy level. This is due to 
the fact that the gas phase is in reality more complex and 
consists of various components at different temperatures. 
The S(4) to S(7) lines (upper energy level from ~3500 
to 7200 K) show in general a much smaller slope than 
the slope inferred from the S(0) and S(l) lines (upper 
energy level of ~500 to 1000 K), representing a "hot" and 
"warm" temperature component of the gas, respectively. 
To take into account these two temperature components, 
we fitted a double exponential profile to our data, given 
as 

— = c + ciexpf Ei/T warm ) + C2exp(Ei/Thot) (1) 
9i 

where iVj is the number of molecules, gi the statistical 
weight of each state, Ei the upper Energy level, C0-2 
constants, and T warm and That the two excitation tem- 
peratures. This fit has the advantage that it simultane- 
ously fits two components, rather than starting with two 
independent temperature estimates using different lines 
and subsequently subtracting the contribution of the hot 
component from the warm component. 

Fig. [TU] shows the excitation diagrams and the fits to 
the warm and hot H2 phases for the different regions 
of our merger systems. For all our objects, except for 
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NGC 6621, the excitation diagram of the entire merger 
system could be well approximated by a two tempera- 
ture component fit. Only the warm component could 
be measured for NGC 6621 because the H2 lines S(4)- 
S(6) are not detected. We find for our sample a mean 
temperature of the warm H2 phase of 255K, but with 
a large variance between different regions, ranging from 
87 K (NGC 520 South) to 348 K (NGC 7592 South). 
The mean temperature of the hot phase for our sample 
is 830 K, ranging from 550K (NGC 520 South) to 1079 K 
(NGC 7592 South). In comparison to other studies, our 
average excitation temperature of the warm component 
seems to be slightly larger than in the nuclear region 
of a sample of SINGS galaxies (~200 K; iRoussel et al.l 
l2007fl . but smaller than in the nuclear region of star- 
burst galaxies (~360 K; Bcrnard-Salas ct al. 2009) and 
ULIRGs (~336 K; iHigdon et al.ll2006f ). Note that this 
comparison does neither take into account the different 
sizes over which the excitation temperature is averaged 
nor the large differences in warm molecular masses in- 
volved. 

We used the two excitation temperatures T warm and 
That to estimate t he to tal warm and hot gas masses 
(Rigo poulou et all l2002t ). M warm and Mhot, hi our 
sources, which are presented in Table |U The total mass 
of the warm component (median M warm = 3.3 x 10 7 M ) 
is in general at least an order of magnitude larger than 
the total mass of the hot component (median M} lot = 
0.1 x 10 7 M ). 

3.3.3. Temperature maps 

To reveal spatial variations of the H2 temperature in 
the merger components, we compute H2 temperature 
maps as shown in Fig. [IT] In only one galaxy, NGC 6240, 
does the H2 line emission from the the hot component 
(S(4)-S(7)) have a sufficient signal-to-noise ratio per pixel 
to employ a two component temperature fit per pixel. 
(Fig. ri2"j) . For all other galaxies a one-component fit (per 
pixel) reveals similar results on the warm component as 
the two-component fit (per pixel) due to a limited num- 
ber of significant H2 lines detected per pixel. Note that 
this does not imply that these galaxies have only one 
temperature component per pixel, since our spatially in- 
tegrated flux measurements (with higher signal-to-noisc) 
show that almost all galaxies have at least two temper- 
ature components (§ I3.3.2|) . Hence, the absolute tem- 
perature scale is not neccessarily equal to the spatially 
integrated temperature, which is due to smaller signal- 
to-noise ratios on pixel-scale. The absolute temperatures 
are given in the previous section (see also Table U]) and we 
focus instead on the temperature distribution in this sec- 
tion. In general, we find an increase in temperature to- 
wards the peak of the combined H2 flux, which is in most 
cases aligned with the center of the main galaxy compo- 
nent. We find no Hi hotspots between merger galaxies. 
An exception is NGC 6240 where we find a faint sec- 
ond temperature peak (500 K), ~7 kpc eastwards, not 
in alignment with the Hi flux intensity nor any other 
mid-IR feature or NIR concentration. Although both 
gas temperature components are mapped, only the "hot" 
component exhibits this secondary temperature peak. 

3.3.4. Excitation Mechanisms: Evidence for an Intrinsic 
Relation between H2 and PAH Emission 



One of the major excitation sources in starburst galax- 
ies is the far-ultraviolet (FUV) radiation of massive stars 
in photo-dissociation region s (PDRs) , with photon ener - 
gies between 6 and 13.6 eV (Hollcnba ch fc Tielensl[T997T ) . 
PDRs are attenuated by the absorption and scattering 
of UV photons from a hot star, generating several layers 
with different temperatures and material components. In 
general, PAH and H 2 can be expected to be co-spatial 
since both spe cies ca n be excited by FUV photons (as 
discussed in § I3.3.1j) . Galactic studies of the molecu- 
lar hydrogen emission in PDRs suggest that the forma- 
tion of H 2 occurs on the surfaces of PAHs (see, e.g., 
lHabart et all 120031: iVelusamv fc Langerl [20081. A simi- 
lar correlation bet ween H2 and PAH emission is seen in 
exter nal galaxies (iRigopoulou et ah [2002; Rousse fet al.l 
120071 ). In particular. IRoussel et al.1 <|2007n found for nor- 
mal star-forming galaxies, on average, a constant flux ra- 
tio of H 2 (S0-S2)/PAH~0.008 over a large range of total 
H 2 fluxes, suggesting relatively constant physical condi- 
tions in PDRs. However, this study was limited a) to 
the average emission of only the H II nuclei and com- 
plexes within the kpc-scale circumnuclear region for each 
galaxy, and b) by the fact that the power emitted in 
the aromatic bands was estimated from the IRAC4 fil- 
ter rather than via direct spectroscopic measurement of 
PAH emission lines. On the other hand, a recent study 
of ULIRGs suggest that the H2 emission does not come 
from PAH regions (PDR's) and originates outside the 
dust obsc ured regions du e to shocks in the surrounding 
material (Zakamska 20"Toh . 

Since H2 emission can arise from at least three differ- 
ent mechanisms (UV-excitation in PDRs, shocks, or hard 
X-ray radiation from an AGN), one important question 
is whether there is one single excitation mechanism that 
dominates the radiation power emitted in the mid-IR or 
a mix of various excitation mechanisms in mergers. In 
principle, studying the relation between PAH and rota- 
tional H2 emission over various spatial scales allows us 
to constrain the role of PDRs, as both, PAH and H2, are 
tightly correlated with the physical conditions in PDRs 
(in particular with the radiation field intensity Go and 
the hydrogen density tih)- Since most of the warm H2, 
at the lowest rotational temperatures, emits primarily in 
the S(0) to S(2) lines, whereas the S(3) line emission has 
a noticeably higher contribution from hotter H2, proba- 
bly indicating more mixed excitation, a useful measure 
of the warm H2 in star-forming galaxie s is the sum of the 
S(0) to S(2) lines (jRoussel et al.ll200l . 

Here we study the relation between the combined H2 
flux (from SO to S2) and the sum of the four main PAH 
bands (at 6.2, 7.7., 11.3, and 17 /im) for various spa- 
tial locations in major merger systems (integrated over 
apertures of 3 x 3 pixels, or 5.6 x 5.6", corresponding 
to ^2 kpc at a mean distance of 79 Mpc). In Fig. [13] 
wc present the total flux ratios H2(S0-S2)]/PAH as a 
function of the total H 2 (S0-S2) luminosity, while Fig.[T4l 
shows the logarithmic relation between H 2 (S0-S2) and 
PAH for each galaxy. The following results are derived, 
and compared to the cen ter of normal star- forming galax- 
ies (jRoussel et al.ll2007i ): 

1 . On global scales (integrated over the entire merger 
component), the majority of galaxies in our 
merger sample have similar total H2(S0-S2)/PAH 
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ratios (mean: 0.008 ± 0.004) as normal star- 
forming galaxies (H 2 (S0-S2)/PAH-0.008). How- 
ever, three galaxies in our sample (NGC 4676 
South, NGC 6240, and NGC 3921) have a signif- 
icant (a factor of ten) larger global ratio (H 2 (S0- 
S2)/PAH=0.08), clearly separated from the rest of 
our sample and normal star-forming galaxies, and 
even larger than most AGN galaxies in the SINGS 
sample (mean H 2 (S0-S2) /PAH=0.016 for Se yfert 
and LINER host galaxies IRoussel et al~ll2007| ). 

2. On local scales (measuring different apertures 
within each merger component), we find an intrin- 
sic relation between H 2 and PAH flux as given by 
a power law, 

log[H 2 (S0 - S2)] = k log[PAH] + log[const.] (2) 

All galaxies in our sample follow this power law 
with a similar slope component k = 0.61 ±0.05 but 
a different constant (up to a factor of 10 difference). 

3. Although the center of galaxies has the strongest 
radiation field and hosts an AGN for at least two 
galaxies in our sample (see Tabled]), it follows the 
same power law as regions in the outer disk. We 
find no clear correlation between the presence of an 
AGN and the central H 2 /PAH ratio. 

The presence of a power law suggests that H 2 (S0- 
S2) and PAH emission are tightly correlated, given as 
H 2 (S0-S2)~ PAH fe , where k is the scaling exponent 
and a a constants. The scale invariance of the power 
law becomes clearer in the logarithmic form, log[H 2 (S0- 
S2)]=fc log[PAH]+log[a], with a mean slope (exponent) 
k = 0.61 (±0.05) for our sample and a mean constant of 
log [a] = —4.0 ± 0.6. Note that the scatter around this 
slope typically increases towards the outer disk (smaller 
signal-to-noise) and that we can not rule out a possi- 
ble dependence of the power law on the resolved spatial 
scales. In a number of galaxies (NGC 6240, NGC 520, 
NGC 2623, and NGC 3921) the faint regions typically fall 
below our best fit power law for the entire galaxy. These 
regions may have a steeper slope, closer to a constant 
value, implying the bright regions may have systemati- 
cally less H 2 emission for a given PAH flux than would be 
extrapolated from a fit to the faint regions in the outer 
disk. A higher extinction for the H 2 lines might explain 
this, but higher signal to noise data is required to confirm 
this change in slope. 

To ensure that the slopes in Fig. [TJ] and the fact that 
the nuclei display lower H 2 /PAH ratios than the galactic 
disks was not the result of having higher spatial resolu- 
tion in the dominant PAH feature (7.7 /xm) compared 
to the dominant H 2 feature (17 /im) because of the 
diffraction-limited nature of Spitzer, we have smoothed 
the SL PAH maps by a factor of three to match SL to 
LL, and re-calculated the H 2 /PAH ratios for all galaxies. 
After this process, the slopes in Fig. IT4l and the ratios in 
the nuclei show a change of less than 3%, suggesting the 
effects we are seeing, namely slopes significantly less than 
unity and lower H 2 /PAH ratios in the nuclei than in the 
disk, are real. Moreover, we find a similar trend for most 
of our galaxies if we calculate the H 2 /PAH ratios using 
the H 2 S(3) line (9.7 fim) instead of H 2 (S0-S2), but with 



a typically much larger scatter and uncertainty (likely 
due to the fact that the S(3) emission line is more af- 
fected by uncertainties in the dust exctinction correction 
and that it has a noticeably higher contribution from 
hotter H 2 with a more mixed excitation as mentioned 
above) . 

How ever, compared to the study of IRoussel et al.l 
(|2007f ). which suggests a linear relation with a constant 
ratio H 2 (S0-S2)/PAH=0.008 for the center of normal 
star-forming galaxies, we find that a constant ratio (i.e. 
k = 1) can not describe the relation between PAH and 
H 2 emission over the entire galaxies' disk as demon- 
strated in Fig. [TB I The nuclei of the galaxies tend to 
be closest to the IRoussel et al.l (|2007l ) value while the 
regions in the disk tend to have higher H 2 /PAH ratios 
(increased H 2 /PAH ratios at radii of 5 to 15 kpc have 
been also found in the edge-on d isk galaxies NGC 4565 
and NGC 5907: lLaine et al.ll20ldh . Instead, the intrinsic 
relation between PAH and H 2 emission follows a power 
law from the center to the outskirts of major mergers. 
This tight correlation and the fact that we find a simi- 
lar slope (exponent) for diverse merger systems suggest 
a common underlying excitation mechanism. The most 
likely mechanism is FUV excitation of PAH and H 2 gas 
in PDRs, since other mechanisms, such as shocks due 
to supernovae or AGN, would a) not show the same be- 
havior over all spatial scales and b) are not expected to 
follow a tight relation such as a power law. The coupling 
between PAH and H 2 implies that both are excited pre- 
dominantly in PDRs, although they may not come from 
the exact same layers (at the same optical depths within 
the clouds). 

A possible explanation for the large range of the con- 
stant a may be caused by the different physical prop- 
erties of the merger systems, e.g. a different abun- 
dance of warm molecular gas to PAHs. This would 
imply that any correlation between H 2 and PAH emis- 
sion based only on global measurements for a sample 
of different galaxies would be washed out, simply be- 
cause the global measurement reflects only the static 
physical properties rather than the behaviour of the H 2 
and PAH excitation within the galaxies. In princible, a 
possible contribution from mechanical heating through 
shocks may play a role as well on some local scales (e.g., 
for the center of NGC 6240, see § [3T1]) . However, this 
is likely not the dominant contribution in our galax- 
ies since observat i ons of extreme H 2 -luminous galaxies 
(lOele et all 120101 : iNesvadba etafl 120101 : iGuillard et all 
120111 ) have shown that shock-excited galaxies have typ- 
ically a larger H 2 /PAH ratic0 than in our sample 
(e.g., the Stephans Quintet with H 2 (S(0)-S(3))/PAH 
7.7 um>l IGuillard et al.llMToh . 

3.4. Local Variations versus Evolutionary Trends along 
the Toomre Sequence 

To investigate possible evolutionary effects on the gas 
and dust properties of major mergers we have studied 
several mid-IR features as a function of their merger 

13 Note that in some studies the H2/PAH ratio is defined as 
H2(S(0)-S(3))/PAH 7.7 (im which is for our sample typically a 
factor of ~3 larger than H 2 (S(0)-S(2))/PAH(6.2, 7.7, 11.3, 17 /im), 
but the former ratio shows typically a larger variation on local 
scales and a less significant power law than the latter one. 
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stage, representing the early (NGC 4676, NGC 7592, 
NGC 6621 - distinguishable parent disks and large 
tidal tails), mid (NGC 2623, NGC 6240, NGC 520 - 
highly overlapping disks and/or double nuclei) and late 
(NGC 3921, NGC 7252 - single nucleus, elliptical-like 
light distribution) stages. For instance, a change in dust 
grain size, ionization stage, neutral to ionized gas frac- 
tion, or gas temperatures could be expected. To test 
such a possible scenario we study the following mid-IR 
diagnostics as a function of merger stage as shown in 
Fig. HU [Ne III]/[Ne II], the PAH interband strength ra- 
tios 11.3/17 /im and 7.7/11 /im, the warm gas to PAH 
ratio H 2 (S0-S2)/PAH as described in § |3"X?I the dust 
continuum ratios 30/15 /im and 15/5.5 /im as derived in 
§ 13.2.11 and the two warm H 2 gas components. The re- 
sults are presented for the entire merger systems as well 
as for the individual merger components. Furthermore 
we explored the radial profile and spatial size (derived 
from Sersic profile fit) of several mid-IR features as a 
function of merger stage, such as the radial size ratio of 
H 2 /IR continuum, PAH 7.7/11.3 /mi, and H 2 /PAH. Our 
main conclusion is that the variety of mid-IR properties 
in different regions as shown in Fig. [16] is significantly 
larger than any trend along the merger sequence using 
the global properties. This does not necessarily imply 
that an overall trend along the merger sequence does not 
exist, for which we find indications in the mid-IR con- 
tinuum slope (Fig. [TBI and the relative radial size of the 
PAH interband emissions and the H 2 emission as shown 
in Fig. 1171 However, a much larger sample would be re- 
quired to confir m statistically such a, trend in the mid-IR 
properties (see iStierwalt et ail 120 111 for nuclear mid-IR 
properties as function of merger stage). 

Despite the small number of galaxies in our sample 
we find two galaxies that are outliers in one or more 
properties: NGC 6240 does not only show a significantly 
larger H 2 /PAH flux ratio (see § I3.3.4[) . but also a larger 
PAH 11.3/17/xm flux ratio as well as a larger ratio of the 
radii H 2 /IR cont, H 2 /PAH and PAH 11.3/7.7/mi than 
most of the other galaxies in our sample. Previous ob- 
servations of the center of this gal axy have revealed two 
nuclei with a separation of 1.8" ([Fried fc Schuld ["19831 ) 
and that the H 2 emission peaks between the two nu- 
clei and is thermally excited vi a slow shocks triggered by 
the collision of the two nuclei (Ivan der Werf et al1ll993t 



iSugai et al.l[l^97l iTecza et al.l 120001 : lArmus et al.H2006[ ) 



This may explain the intense H 2 emission in the very 
center of NGC 62 40 (unresolved with Spitzer IRS, see 
lArmus et al.ll200"6T) . but is likely not responsible for the 
correlation seen in the outer disk between H 2 and PAH. 
Another galaxy, NGC 2623, is significantly larger in the 
flux continuum ratio (30/15/xm= 16.3) than the rest of 
the sample (mean 30/15 /im~ 6), which may be due to 
the star-formation intensity peak that occurs during a 
relatively short timescale (< 0.5 Gyrs) at the final coa- 
lescence of t he nuclei as shown in recent major merger 
simulations ([Cox et all 120081; iHopkins et all 120091) . Fur- 
thermore, this stage also seems to show a relative change 
in the radii of the mid-IR features, namely, a maximum 
in the H 2 emission line radius to mid-IR continuum ra- 
dius (at 5.5-15 /jm) as well as in the PAH 11.3 /um to 
7.7 /im line ratio radius. 



3.5. Multiwavelength View of the Interstellar Material 
in Early Merger Stages 

In this section we discuss in more detail three "clas- 
sical" major mergers, NGC 4676, NGC 7592, and 
NGC 6621, and compare their mid-IR properties with the 
distribution of the stellar (NIR a nd FUV), cold atomic 
(HI), and cold molecular (CO, see lYun fc Hibbardl[200l 
gas components. These three mergers are in particu- 
lar interesting since they are caught in the early merger 
stages and have spatially resolved components which al- 
lows us to compare in detail the physical properties be- 
tween the two nuclei in each merger system. 

3.5.1. NGC 4676 

Nicknamed the "Mice", each of the two galaxies in 
NGC 4676 possesses a long tail as observed in the op- 
tical light. The two galaxies of this merger system have 
a projected separation distance of 36.8"(~17 kpc). In 
particular, the northern galaxy exhibits a long tidal tail 
with a projected extent of 48 kpc (measured in the NIR) , 
visible in the FUV, near-infrared, mid-IR and HI emis- 
sion (see Fig. [T5[) . Material is being exchanged between 
the two disks as indicated by the pre s ence o f a b right op- 
tical b ridge between them. iStocktonl ([1974D and lHibbardl 
(1995) found that the tails have widespread star forma- 
tion which accounts for 16% of the total Ha emission 
and exhibit a high atomic gas content with on average 
Mhi/Lr ~ 0.4 MqLq 1 , with the stellar tails exceeding 
the gaseous tails in extent. 

Although both galaxies have roughly the same near-IR 
luminosities (as measured in the H-band), we find signifi- 
cant differences in the distribution of the mid-IR features: 
As shown in Fig. [TH] most of the PAH and CO emission 
(each >90%) originate from the northern galaxy, while 
the main H 2 line emission (>80%) is found in the south- 
ern galaxy. This finding may suggest a significantly dif- 
ferent stellar population or distribution of molecular gas 
and young stars in the merger components. The north- 
ern component has a significantly smaller [Ne III]/[Ne 
II] ratio (0.13) and larger PAH 7.7/11.3 fim ratio (4.5) 
than the southern galaxy (~0.4 and 3.5, respectively), 
suggesting that the average stellar population is younger 
and the grains may be more inonized in the northern 
component. Interestingly, we find that the H 2 tempera- 
ture components (^300 K for the warm and ^1000 K for 
the hot component) are roughly the same in both galax- 
ies. Thus, the strong H 2 emission in the southern galaxy 
may be simply due to the presence of a larger amount of 
warm molecular gas (H 2 mass two times larger than in 
the northern galaxy). 

Furthermore, we find PAH emission in the northern 
tidal tail, but no significant H 2 emission. Fig. [19] shows 
that the interband strength of the PAH 7.7/11 /mi and 
6.2/7.7 /im ratios of the tidal tail regions are all char- 
acterized by the ionized phase (using the model by 
iDraine fc Lill200ll ). This is perhaps not surprising given 
the large number of young stars in the tidal tail as seen 
in the strong FUV emission (F ig. [TBI and Ha emission 
(IStocktonl [1971: iHibbardl IT995I: iHibbard fc van Gorkoml 
1996)), demonstrating that the grain properties in inter- 
acting galaxies as traced by the mid-IR feature ratios can 
be affected over large scales in interacting systems. 
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3.5.2. NGC 7592 

The two mergi ng; galaxies in thi s LIRG system 
(Lir/Lq = 10 1133 , Sa nders et al.ll2003f ) have a projected 
separation of 14" (~7 kpc), a NIR luminosity ratio of 
-1:1, and show tidal tails in the NIR light. Fig. [20] high- 
lights the details of this system at multiple wavelengths. 
Like NGC 4676, the PAH and H2 emissions have a differ- 
ent distribution, with the PAH emission primarily origi- 
nating from the eastern starbursting component and the 
H2 from the western comp onent, which has a Seyfert 2 
nucleus (jLaine et al.l 120031 ). 

We find an additional concentration in the overlap re- 
gion between the two nuclei (with an offset of —12" 
towards south) in the NIR light, FUV light, HI emis- 
sion and in several mid-IR features (main PAH fea- 
tures, [Nell] and [Nelll]). No significant H 2 emission 
is found at this position. The PAH interband strength 
and [Nelll] / [Nell] ratios indicate that most of the emis- 
sion from this intermediate region originates from ionized 
gas (see Fig. [7]), which is likely due to star-formation 
as we find a stron g emission peak in the FUV and Ha 
(jLaine et al.l l2003]) at this position as well. 

3.5.3. NGC 6621 

The two merging components of NGC 6621 have a pro- 
jected separation of 41"(— 17 kpc). Unlike the two previ- 
ous discussed major mergers, NGC 4676 and NGC 7592, 
which have a NIR light ratio of 1:1, the two components 
of this system have a ratio of —1:2. However, this ratio 
changes dramatically in the light of different wavelengths. 
The mid-IR light of the smaller component, as given by 
our IRS continuum measurement, is at least two orders 
smaller than the larger component (nearly all of the mid- 
IR emission in NGC 6621 is coming from only one of the 
two nuclei). As shown in Fig. I21[ we find no significant 
FUV, H2, CO, and HI emission at the position of the 
smaller galaxy, and only very faint emission from PAH 
features. Instead, we find significant emission originating 
from a peak in the overlap region between the two merg- 
ing galaxies. The interband PAH strength ratios of this 
peak (see Fig. [7]) indicate that the PAH component may 
originate partially from neutral gas layers. However, we 
also find significant H2 and FUV emission in this overlap 
region, suggesting a burst of star-formation here. 

The distribution of the c old molecular gas , as observed 
in the CO emission line (|Iono et al.l 120051) , is not co- 
spatial to the stellar (NIR), warm gas (H2) and PAH 
distribution: besides the CO emission peak of the main 
galaxy, we find an additional CO concentration (—30% 
of the total CO emission) between the main galaxy and 
the overlap region as seen in the FUV and NIR light, 
and mid-IR features. This significant offset suggests a 
large reservoir of cold molecular gas with no signs of star 
formation in the mid-IR features. 

4. CONCLUSIONS 

We analyzed the spatial distribution and properties 
of the main mid-IR features observed with Spitzer/IRS 
for eight major merger systems of the Toomre Sequence. 
Maps of the key diagnostic features of the physical condi- 
tions in these interacting starburst galaxies are presented 
— such as PAH ratios, ionized emission line regions, the 
warm molecular gas (H2) distribution and temperatures. 



Moreover, we have compared the mid-IR features to the 
cold atomic and molecular gas distribution to obtain a 
comprehensive view about the variation and evolution of 
the interstellar material phase space in major mergers. 
The main results are summarized in the following: 

1. The range and distribution of the global PAH in- 
terband strength and [Ne III] / [Ne II] ratios of ma- 
jor mergers are very similar to normal star-forming 
galaxies. However, regarding the distribution of 
these properties over all spatial scales within the 
galaxies, we find significantly different PAH popu- 
lations (defined by their ionization state and grain 
size) among the merger systems in our sample. 

2. The cold molecular gas component as observed in 
the CO emission line does not usually correlate 
with the warm gas in early merger systems, leading 
to large spatial variations in the ratio of cold-to- 
warm molecular gas. In particular, the cold molec- 
ular gas distribution in NGC 6621 shows a signifi- 
cant offset not only from the warm gas component, 
but also from the stellar (NIR) and dust distribu- 
tion (PAH, mid-IR continuum). 

3. Two merger systems (NGC 7592 and NGC 6621) 
show evidence for star formation in overlap regions 
between the nuclei, namely strong PAH, [Ne II] , [Ne 
III] and warm gas (H2) emission, but no apparent 
corresponding CO (J=l-0) emission. 

4. The majority of galaxies in our merging sequence 
sample have similar global H2-to-PAH flux ratios 
as normal star-forming galaxies with H2 temper- 
atures ranging from — 90K to HOOK. Three star- 
burst galaxies in our sample (NGC 4676 South, 
NGC 6240, and NGC 3921) have a factor of ten 
larger ratios, clearly separated from the rest of our 
sample and normal star-forming galaxies, and even 
larger than wha t it is in most AGN galaxies in the 
SINGS sample (|Roussel et al.ll2007T ). However, re- 
garding the different apertures (locations) within 
a merger, we find an intrinsic relation between H2 
and PAH flux as given by a power law with a very 
similar exponent (0.61 ± 0.05) for all galaxies but 
with a constant that can vary up to a factor of ten 
between different galaxies. Such a tight correla- 
tion between PAH and H2 over all spatial locations 
suggests a common dominant excitation mecha- 
nism for H2 emission in major mergers, which is 
most likely absorption of UV photons in photo- 
dissociation regions rather than shock excitation, 
since PAHs are known to arise predominantly from 
these regions. 

5. The minimum of the 11.3/7.7 /Ltm PAH interband 
strength ratio is typically located in the nuclei of 
galaxies, suggesting a larger ionization state of the 
ISM towards the center, as expected from cen- 
tral starbursts or low-luminosity AGNs (Seyfert, 
LINER). The [Ne III/[Ne II] ratio increases with 
distance from the nucleus, seemingly conflicting the 
expectation of a decrease of the hardening of the 
radiation field with distance from the center. A 
possible explanation for a suppressed [Ne III]/[Ne 
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II] ratio in the center is either a larger contribu- 
tion of the old stellar population in the nuclei, or 
a fraction of the most massive stars are hidden in 
ultra-compact HII regions due to high dust densi- 
ties in the nuclear regions. Thus, the [Ne III/[Ne II] 
ratio is likely more sensitive to the average stellar 
age and/or gas density rather than a hardening of 
the field, suggesting that the 11.3/7.7/um PAH ratio 
is a better diagnostic of the ionization state in all 
merger components, as seen in the center and out- 
skirts of galaxies, in overlapping regions between 
merger components, and in tidal tails. 

6. Variations of the internal (local) physical condi- 
tions within a merger are much larger than the 
systematic (global) variations along the Toomre Se- 
quence. Note that this does not necessarily exclude 
the possibility of an evolutionary trend, for which 
we find indications in the mid-IR continuum slope 
and the relative radial size of the PAH interband 
emissions and the H 2 emission. 
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stitute of Technology, under contract with the National 
Aeronautics and Space Administration. We thank the 
referee for a careful review and valued suggestions, and 
Pierre Guillard for a useful discussion. 
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(7) NGC3921 (8) NGC7252 
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Fig. 1. — Spitzcr IRAC 3.6 /im images of our sample along the Toomre Sequence (from 1-8, top-left to bottom-right). The Field-of-View 
of the Spitzer IRS SL images is overlaid in rectangles and the figures are oriented such that north is up and east left. 
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FlG. 2. — Mid-IR continuum maps of our sample (see § 13. II for more details). The dashed-line rectangles show the individual regions for 
which spectra have been extracted. North and east are marked by arrows in each frame and the bar corresponds to 5 kpc. 
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Fig. El— (Continued). 
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Fig. 4. — The main mid-IR feature maps with contours, from left to right and top to bottom: the integrated mid-IR continuum (5.3- 
15 fim) as measured with PAHFIT, the combined PAH features, the combined H2 emission features, and the [Ne II] emission (see text). 
The center is marked with a plus symbol signifying the peak in the integrated mid-IR continuum and the tickmarks correspond to 1 
pixel-intervals (1.85"). The North and East directions in each frame of this figure are the same as in the corresponding frame of Fig. [2] 
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Fig. [4] — Continued. 
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Fig. 5. — Ratio map of PAH 6.2/7.7 (im overlaid with the combined PAH emission (6.2,7.7 and 11.3 /an) in contours. The North and 
East directions in each frame of this figure are the same as in the corresponding frame of Fig. [2] 
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Fig. [7] — Continued 
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Fig. 8. — [Nc III]/[Nc II] line flux ratio maps for the Toomre sequence galaxies. The contours indicate the sum of the [Nc III] and 
II] flux. The North and East directions in each frame of this figure are the same as in the corresponding frame of Fig. [2] 
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Fig. 9. — H2 feature maps with contours (from S(5) at top left to S(0) at bottom right). The plus marker indicates the central peak in 
the mid-IR continuum emission. Figures for the rest of the sample are available in the online version of the Journal. The North and East 
directions in each frame of this figure are the same as in the corresponding frame of Fig. [2] 
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Fig. 10. — H2 excitation diagram for NGC 6240, showing the upper level population divided by the level degeneracy as a function of 
upper level energy (for the lines observed with the IRS). The lines are measured with PAHFIT, including a correction for extinction. 
Uncertainties are indicated as vertical bars for each line. A double exponential profile fit (solid line) takes into account two temperature 
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Fig. 11. — The single-component H2 temperature map (in colors) with the combined H2 emission (in contours) for NGC 4676. Figures 
for the rest of the sample are available in the online version of the Journal. 
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Fig. 12. — The two H2 temperature maps (left: hot, right: warm component) for NGC 6240, with the combined H2 emission in contours. 
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Fig. 13. — Ratio of the power emitted in the sum of the H2 S(0) to S(2) transitions to the power emitted in the main PAH bands (at 
6.2, 7.7, 11.3, and 17 /jm) over a large range of H2 luminosities. Measurements of the entire merger system are marked with filled squares 
while the individual merger components are marked with open circles. The solid (dotted) and dashed (dashed-dotted) line represents the 
mean (uncertainty) H2( S0— S2 ) /PAH value measured for the center of normal star-forming galaxies and LINER/Seyfert nuclei in the SINGS 
sample, respectively. ( Rousscl ct al. 2 0071 ) . 
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Fig. 14. — The power emitted in the sum of the H2 S(0) to S(2) transitions as a function of the power emitted in the main PAH bands 
(at 6.2, 7.7, 11.3, and 17 (im) in logarithmic scales and in units of W m _2 sr _1 . The crossed circles represent the central aperture of a 
merger component. The distribution is fitted using a power law (solid line), given as log[H2(S0— S2]= fcxlog[PAH]+c with the slope k 
and the constant c. The dashed line represents the constant ratio model for the center of normal star-forming galaxies, given as H2(S0- 
S2)/PAH=0.008 URoussel et al.|[20Tj7f) . 
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Fig. 15. — Overview of the power laws (log[H2(S0-S2]= fcxlog[PAH]+c) as derived in Fig |14l for each galaxy (solid li nes) . The dashed lin e 
represents the constant ratio model for the center of normal star-forming galaxies, given as H2(S0-S2)/PAH=0.008 (Rousscl ct al. 2007), 
clearly demonstrating that a constant ratio (i.e. k = 1) can not describe the measured relation between PAH and H2 (mean k = 0.61±0.05). 
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Fig. 16. — Dust and gas properties as a function of the Toomrc Sequence (from the early merger stage NGC 4676 to the late stage 
NGC 7252). The entire merger systems and individual merger components are marked as filled squares and open circles, respectively. From 
left to right and top to bottom: ratio [Ne III]/[Ne II], ratio of the sum of the H2 S(0) to S(2) transitions to the main PAH bands, PAH 
ratio 11.3/7.7 fim, PAH ratio 7.7/11.3 fim, the mid-IR continuum ratio 30/15 fim, 15/5.5 fim, the "warm" H2 temperature, and "hot" H2 
temperature component (see text for details). 
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Fig. 17. — The relative radii of the major merger components as a function of the Toomre Sequence (from the early merger stage 
NGC 4676 to the late stage NGC 7252). Left: The ratio of the strongest H2 emission line (S(3)) to the mid-IR continuum radius (at 
5.5—15 /im) as a function of merger stage. Right: The ratio of the PAH 7.7 /an to 11.3 )im line radius as a function of merger stage. 




Fig. 18. — Multiwavclength view of NGC 4676, from left to right and top to bottom: The combined main PAH features, the warm 
molecular gas (H2 emission), the cold atomic gas (HI), the cold molecular gas (CO), the NIR (Spitzer IRAC 3.6 (im), and the FUV 
(GALEX). The North and East directions in each frame of this figure are the same as in the corresponding frame of Fig. [2] 
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Fig. 19. — Variation of interband strength ratios of the PAH features in the tidal tail of NGC 4676 indicates emission from ionized PAHs. 
Same scales, markers and labels as in Fig. 
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Fig. 20. — Multiwavelength view of NGC 7592, from left to right and top to bottom: The combined PAH features, the warm molecular 
gas (H2 emission), the cold atomic gas (HI), the cold molecular gas (CO), the NIR (Spitzer IRAC 3.6 /^m), and the FUV (GALEX). The 
North and East directions in each frame of this figure are the same as in the corresponding frame of Fig. [2] 
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Fig. 21. — Multiwavclength view of NGC6621, from left to right and top to bottom: The combined PAH features, the warm molecular 
gas (H2 emission), the cold atomic gas (HI), the cold molecular gas (CO), the NIR (Spitzer IRAC 3.6 /im), and the FUV (GALEX). The 
centers of the two galaxies (as seen in the NIR) and the star-formation peak in between are indicated with a blue plus marker. The North 
and East directions in each frame of this figure are the same as in the corresponding frame of Fig. [2] 
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TABLE 1 

Sample Overview 



Name 


RA 


DEC 


D 




T.S. 


AGN type 




(J2000) 


(J2000) 


[Mpc] 


logio [Lq] 






NGC 520 


01:24:35.1 


+03:47:33 


27.2 


11.91 1 


5 


HII 


NGC 2623 


08:38:24.1 


+25:45:17 


80.6 


11.54 1 


6 


LINER Sy2 


NGC 3921 


11:51:06.9 


+55:04:43 


84.5 


<10.35 2 


7 


HII 


NGC 4676 


12:46:10.7 


+30:43:38 


96.0 


10. 88 3 


1 




NGC 6240 


16:52:58.9 


+02:24:03 


103.0 


11. 85 1 


4 


LINER Sy2 


NGC 6621 


18:12:55.3 


+68:21:48 


85.2 


11. 23 1 


3 


HII 


NGC 7252 


22:20:44.7 


-24:40:42 


62.2 


10.75 3 


8 




NGC 7592 


23:18:22.2 


-04:25:01 


97.2 


11. 33 1 


2 





Note. — Summary of the properties of our sample of merger systems. Column 
(1): Source Name from NASA IPAC Extragalactic Database (NED), Column (2): 
right ascension from NED (J2000), Column (3): declination from NED (J2000), 
Column (4): The luminosity distance in Mpc (adopting Ho = 70 km s — 1 Mpc), as 
pro vided by NED , Colu mn (5): The total infra red l uminosity in l ogio Solar units 
(1: ISanders et all | |2003T) , 2: IChitre & j oe (2002), 3: Brandl et al. (200J)), Column 
(6): The merger stage of the Toome Sequence (from 1—8), Column (7): The AGN 
type derived from the optical and UV spectra, as recorded in NED. 



TABLE 2 

Mid-IR Emission Line and PAH Features 





Name 


RA 


DEC 


Reg 


ion Size 


PAH 6.2/mi 


PAH 7.7/xm 


PAH 11.3itm 


PAH 17itrn 


[Ne II] 12.8lAtm 






(J2000) 


(J2000) 


10" 
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NGC 


520 System 


01:24:35.10 


+03:47:33.0 




75.5 


389.5 ± 2.4 


1427.7 ± 7.1 
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34.5 ± 0.2 


NGC 


520 North 


01:24:33.36 


+03:48:02.5 
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± 4.8 


264.5 ± 14.2 
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520 South 


01:24:34.90 


+03:47:29.5 
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± 
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2623 
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33.5 
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4676 System 
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830.9 ± 7.7 
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4676 South 
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6621 South 


18:12:59.78 


+68:21:14.8 




51.5 






38.0 ± 1.8 


19.9 


± 


0.5 


9.0 


± 


1.1 


0.8 ± 0.1 


NGC 


6621 Middle 


18:12:58.58 


+68:21:28.9 




10.5 


135.1 


± 7.0 


142.9 ± 3.4 


55.9 


± 


1.8 


34.3 


± 


2.3 


4.5 ± 0.2 


NGC 


7252 


22:20:44.77 


-24:40:41.7 




34.2 


188.6 


± 2.6 


660.6 ± 8.8 


191.8 


± 


3.1 


92.5 


± 


2.9 


18.4 ± 0.2 


NGC 


7592 System 


23:18:22.20 


-04:25:01.0 




27.6 


321.5 


± 3.4 


1146.4 ± 13.0 


250.9 


± 


1.9 


127.3 


± 


2.0 


35.8 ± 0.2 


NGC 


7592 West 


23:18:21.84 


-04:24:57.1 




5.3 


350.8 


± 6.5 


1422.0 ± 14.9 


297.4 


± 


2.3 


171.1 


± 


3.4 


53.3 ± 0.4 


NGC 


7592 East 


23:18:22.64 


-04:24:58.1 




5.3 


726.6 


± 6.8 


2446.7 ± 18.6 


508.8 


± 


4.6 


224.3 


± 


4.1 


67.2 ± 0.6 


NGC 


7592 South 


23:18:22.16 


-04:25:08.6 




1.1 


36.5 


± 2.5 


132.3 ± 6.2 


29.5 


± 


1.3 


23.5 


± 


4.0 


4.2 ± 0.2 



Note. — Overview of the integrated PAH and Mid-IR feature strengths for all merger components. Units are in 10" 9 W m 2 sr 



Spitzer IRS Spectral Mapping of Major Mergers 



TABLE 3 
H 2 Fluxes 





Name 


H 2 S(0) 


H 2 S(l) 


H 2 S(2) 


H 2 S(3) 


H 2 S(4) 


H 2 S(5) 


H 2 S(6) 


H 2 S(7) 






28.22^tm 


17.03/im 


12.28/im 


9.66/im 


8.03/im 


6.91/mi 


6.11/im 


5.51/xm 


NGC 


520 System 


5.4 ± 


0.3 


11.2 


± 


0.2 


3.4 


± 


0.2 


10.9 


± 


0.7 


4.9 


± 


0.5 






2.8 ± 0.7 


NGC 


520 North 


0.7 ± 


0.1 


1.9 


± 


0.3 


0.7 


± 


0.4 


0.7 


± 


0.4 


2.6 


± 


1.2 






1.8 ± 1.7 


NGC 


520 South 


96.5 ± 


4.1 


86.4 


± 


0.8 


36.8 


± 


0.6 


139.4 


± 


3. 1 


35.3 


± 


2.0 


21.5 ± 4.1 




3.3 ± 2.9 


NGC 


2623 


1.5 ± 


0.3 


5.5 


± 


0.2 


1.1 


± 


0.2 


6.5 


± 


0.8 


1.9 


± 


0.6 






3.2 ± 0.9 


NGC 


3921 


0.7 ± 


0.1 


5.1 


± 


0.3 


1.6 


± 


0.1 


3.5 


± 


0.5 


0.7 


± 


0.2 


1.4 ± 0.4 


0.5 ± 0.4 


0.2 ± 0.4 


NGC 


4676 System 


0.5 ± 


0.1 


2.9 


± 


0.3 


1.1 


± 


0.1 


1.7 


± 


0.2 


0.8 


± 


0.3 


1.2 ± 0.4 




0.6 ± 0.4 


NGC 


4676 North 


0.3 ± 


0.2 


5.1 


± 


0.6 


2.5 


± 


0.2 


4.1 


± 


0.4 


2.6 


± 


0.5 


2.9 ± 1.0 






NGC 


4676 South 


1.6 ± 


0.2 


9.6 


± 


0.5 


3.6 


± 


0.3 


6.7 


± 


0.4 


1.6 


± 


0.8 


4.3 ± 1.2 




2.2 ± 1.0 


NGC 


6240 


16.9 ± 


0.8 


87.7 


± 


0.3 


25.1 


± 


0.2 


103.3 


± 


0.6 


36.0 


± 


0.6 


63.8 ± 0.8 


12.3 ± 0.7 


28.4 ± 0.8 


NGC 


6621 System 


0.5 ± 


0.0 


2.0 


± 


0.1 


0.9 


± 


0.1 


1.2 


± 


0.1 












0.9 ± 0.4 


NGC 


6621 North 


2.0 ± 


0.1 


5.8 


± 


0.3 


8.8 


± 


0.3 


7.6 


± 


0.4 


6.7 


± 


0.9 


17.2 ± 1.4 






NGC 


6621 South 


0.9 ± 


0.1 


3.5 


± 


0.3 


0.3 


± 


0.2 


2.7 


± 


0.2 


1.3 


± 


0.5 


0.6 ± 0.2 






NGC 


6621 Mid 


0.7 ± 


0.2 


2.3 


± 


0.4 


0.6 


± 


0.2 


0.4 


± 


0.3 












5.0 ± 2.4 


NGC 


7252 


0.4 ± 


0.2 


4.8 


± 


0.3 


2.0 


± 


0.2 


3.7 


± 


0.4 


1.0 


± 


0.6 


0.9 ± 0.8 






NGC 


7592 System 


1.8 ± 


0.1 


7.4 


± 


0.2 


3.8 


± 


0.2 


7.4 


± 


0.5 


1.8 


± 


0.8 


2.6 ± 1.1 






NGC 


7592 West 


2.4 ± 


0.2 


14.1 


± 


O.i 


9.1 


± 


0.4 


13.0 


± 


1.0 


6.8 


± 


1.3 


9.6 ± 2.0 




1.3 ± 1.6 


NGC 


7592 East 


1.6 ± 


0.2 


7.3 


± 


0.1 


1.3 


± 


0.5 


11.5 


± 


1.1 








6.5 ± 2.7 






NGC 


7592 South 


1.8 ± 


0.3 


2.1 


± 


0.7 


0.2 


± 


0.2 


0.2 


± 


0.4 


0.4 


± 


0.7 


0.7 







Note. — Overview of the integrated H2 line strengths of the H2 components for all merger components. The center and size of the 
regions over which is integrated arc the same as in Tab. [2] All units are in 10 9 W m 2 sr 1 . 



TABLE 4 
H 2 Temperatures and Masses 



Name 


T warm 


Thot 


H 2 Mass warm 


H 2 Mass hot 




[K] 


[K] 


[10 6 M ] 


[10 6 M ] 


NGC 520 System 


221. 4± 10.5 


877. 1± 42.3 


16. 0± 0.2 


1.2± 0.3 


NGC 520 North 


246. 1± 93.7 




0.2± 0.1 




NGC 520 South 


87.7± 0.2 


550.1± 7.0 




10.3± 8.9 


NGC 2623 


162. 6± 38.8 


969. 7± 67.8 


118.3± 5.1 


2.6± 0.7 


NGC 3921 


241. 7± 17.4 


664.5± 82.5 


33.4± 1.9 


5.1± 10.0 


NGC 4676 System 


286. 5± 18.7 


1041.8± 188.0 


21. 0± 1.9 


1.0± 0.6 


NGC 4676 North 


309. 4± 27.5 


1091.8± 366.2 


4.9± 0.6 




NGC 4676 South 


291. 7± 17.0 


990.2± 136.7 


9.2± 0.5 


0.7± 0.3 


NGC 6240 


184.9± 1.9 


746.1± 3.2 


1608.1± 6.2 


279. 5± 8.0 


NGC 6621 System 


316. 8± 16.8 




43.8± 1.8 




NGC 6621 North 


380.8± 6.0 




8.6± 0.4 




NGC 6621 Mid 


252. 8± 27.6 




2.9± 0.5 




NGC 7252 


296. 1± 37.0 


663.7± 240.8 


6.5± 0.4 




NGC 7592 System 


333. 8± 26.2 


847.3± 318.4 


14.9± 0.5 




NGC 7592 West 


348. 1± 14.0 


1078.9± 213.3 


5.0± 0.1 


0.1± 0.2 


NGC 7592 East 


339. 4± 66.3 


838.2± 305.8 


2.7± 0.2 





Note. — 



Temperatures and masses of the H2 emission derived by fitting the excitation diagram with a two-component model. 



